Buckled nano rod - a two state system and its dynamics using system plus reservoir 

model 



Aniruddha Chakraborty, 1,2 
department of Inorganic and Physical Chemistry, 
Indian Institute of Science, Bangalore, 560012, India, 
2 Jawaharlal Nehru Center for Advanced Scientific Research, Bangalore, 560064, India. 

(Dated: December 10, 2009) 

We consider a suspended elastic rod under longitudinal compression. The compression can be used 
to adjust potential energy for transverse displacements from harmonic to double well regime. As 
compressional strain is increased to the buckling instability, the frequency of fundamental vibrational 
mode drops continuously to zero (first buckling instability). As one tunes the separation between 
ends of a rod, the system remains stable beyond the instability and develops a double well potential 
for transverse motion. The two minima in potential energy curve describe two possible buckled states 
at a particular strain. From one buckled state it can go over to the other by thermal fluctuations 
or quantum tunnelling. Using a continuum approach and transition state theory (TST) one can 
calculate the rate of conversion from one state to other. Saddle point for the change from one 
state to other is the straight rod configuration. The rate, however, diverges at the second buckling 
instability. At this point, the straight rod configuration, which was a saddle till then, becomes 
hill top and two new saddles are generated. The new saddles have bent configurations and as rod 
goes through further instabilities, they remain stable and the rate calculated according to harmonic 
approximation around saddle point remains finite. In our earlier paper classical rate calculation 
including friction has been carried out [J. Comput. Theor. Nanosci. 4 (2007) 1], by assuming that 
each segment of the rod is coupled to its own collection of harmonic oscillators - our rate expression 
is well behaved through the second buckling instability. In this paper we have extended our method 
to calculate quantum rate using the same system plus reservoir model. We find that friction lowers 
the rate of conversion. 



INTRODUCTION 

Considerable attention has recently been paid to two-state nano- mechanical systems P, 0, d, 0, H, El and the possibility 
of observing quantum effects in them. In the experiments of Rueckes et al. Q crossed carbon nano-tubes were 
suspended between supports and the suspended element was electrostatically flexed between two states. Roukes et al. 
[H propose to use an electrostatically flexed cantilever to explore the possibility of macroscopic quantum tunnelling 
in a nano-mechanical system. Carr et al. [B|, 0] suggest using the two buckled states of a nanorod as the two states 
and investigate the possibility of observing quantum effects. Here we consider a suspended elastic rod of rectangular 
cross section under longitudinal compression. The compression is used to adjust the potential energy for transverse 
displacements from the harmonic to the double well regime as shown in the Fig[T] As the compressional strain is 
increased to the buckling instability the frequency of the fundamental vibrational mode drops continuously to 
zero. Beyond the instability, the system has a double well potential for the transverse motion. The two minima in the 
potential energy curve describe the two possible buckled states at a particular strain 0] and the system can change 
from one state to the other, under thermal fluctuations or quantum tunneling. In our earlier publications 13, II 1 



we have used multidimensional transition state theory (TST) to derive expressions for the transition rate from one 
potential well to the other. We now include the effect of friction on the reaction rate. For this, we follow the procedure 
[12j | used for the study of barrier crossing and other dynamical problems in the presence of friction. 



THE MODEL 



We consider the normal modes and associated quantum properties of an elastic rectangular rod of length L, width 
w and thickness d (satisfying L » w » d) @, El, 14, 15 1. We assume that d is smaller than w so that transverse 
displacements y(x) only occur in the "d " direction. F is the linear modulus (energy per unit length) of the rod and is 
related to the elastic modulus Q of the material by F — Q w d. The bending moment k is given by n 2 = for a rod 
of rectangular cross section and fj, = is the mass per unit length. We take the length of the rod (uncompressed) to 
be Lq. As in the Fig. [2J we apply compression on the two ends, reducing the separation between the two to L. Then 
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FIG. 1: Potential energy V as a function of the fundamental mode displacement Y. The shape of the potential energy is 
harmonic for e > e C) quartic for e — e c = critical strain (e e < 0) and a double well for e < e e . 



if y(x) denotes the displacement of the rod in the '<£' direction, the total potential energy is given by loL[ll| 

L L 

V[y(x)} = \J dx{Fn\y"f + Fe{y'f) + ^{J dx{y'f) 2 + ^V. (1) 
o o 

In the above e = ^ L ~^°- ) is the strain, negative if compressive. We imagine that each segment of the rod (labelled by 
x) is coupled to a collection of harmonic oscillators which form the reservoir. Each section has its own independent 
collection of harmonic oscillators. So the total system, rod plus the reservoir has the energy [10 ] 



N 

H = ^Jd Xf iy 2 +V[y] + \j2jdx 



x m a (x) £ +m a (x)ul(x) ~ ^^K ) • (2) 

[ V m a (x)ujl(x)J J 

In the above, £, a (x) denotes the position of (a, x)-th harmonic oscillator of the reservoir, coupled to y(x), the dis- 
placement of the rod at location x. It has a frequency uj a (x) and mass m a (x), c a (x) determines the coupling of the 
£ Q (x) to the y(x). The way, the coupling has been chosen, the barrier height remains unchanged. Here a varies over 
the collection of harmonic oscillators, and we let a vary from 1 to JV. 
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FIG. 2: The rod under compression: The central figure (A) shows the uncompressed rod of length Lq. On compressing to 
length L, the rod buckle, either to B- or to B+. 

QUANTUM TRANSITION STATE THEORY IN PRESENCE OF FRICTION 

To derive the quantum rate we use the following methodology. The Hamiltonian given in Eq. @ may be treated as 
a quantum Hamiltonian. For a finite discrete set of oscillators one may evaluate the quantum rate using quantum 
transition state theory under harmonic approximations. The quantum transition state theory rate expression is given 
bySEll 
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fric 

quantum 



kT 



h V2sin(M/3Q/ r 



Q* 



-0E a 



(3) 



In the above Q* and Q are the partition function at the saddle point and at the reactant, and iV" c is the frequency 
of the 'unstable mode'. Following Pollak [lil j. the quantum partition function at the saddle point 



N 

n 



L 2 l 2siim(|^xn) fJi 2sinh(±?i/3A* jJ .) 



(4) 



At the saddle point we have N x N — 1 real oscillators and one imaginary oscillator. A„j is the 'bath' frequency and 
Xn is the 'system' frequency at the saddle point. Note the well known divergence of at low temperatures. The 
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quantum partition function at the reactant geometry 

N 



Q 2 ( 2sinh(i?i/3 X „) A 2sinh(±fi/?A nj -) I ' (5) 



In the above A nj - is the 'bath' frequency and Xn is the 'system' frequency at the reactant geometry. 

THE EXTREMA OF THE FUNCTIONAL V[y(x)] 

1y{x 



To find the equilibrium state, we extremise potential energy functional with respect to y(x). For this we put sv ^ y ( x ^ — 



0. This leads to the differential equation 

o 

and the hinged end points have boundary conditions y(0) = y(L) = and y"(0) = y"{L) = 0. The only solution to 
Eq. ([6]) is y(x) — if e > e c , where e c = — . But if e < e c two more solutions are possible. They are given by 

y(x) = ±Ay^sin(^x). (7) 
On substituting Eq. ([7]) into Eq. ©, we find a nonlinear equation for A, which on solution gives 



A = f-^(e c -s). (8) 

These two are the buckled states which are minima of the potential energy surface for e < e c , as all the normal modes 
around this are stable. The solution y(x) = is now a saddle point (see next section for details). 

NORMAL MODES OF THE SYSTEM HAMILTONIAN 

We find the potential energy for the buckled states Ey, = — -^^(e — e c ) 2 . y[x) = is the straight rod configuration 
and is the saddle point for transition from one buckled state to the other. Its potential energy E^J^l is zero. One 
can calculate the barrier height for the process of going from one buckled state to the other as 

&E L B ZZr = ^{e-e c ) 2 . (9) 

The kinetic energy of the rod is ^ y 2 dx, where fj, = is the mass per unit length. We now look at small amplitude 
vibrations around each extremum y ex t(x), which could be the buckled or the saddle point. One gets 

d 2 8y(x,t) } J ( 6 2 V[y(x)} 

On evaluation we get 



= -J dX1 {syJ^SyV)) yM Syix1 ^ (10) 







d 2 Sy(x,t) 

M Op = -CextOy{x,t), (11) 

where the operator C ex t is defined by 

L 

C ext 8y(x,t) =F K 2 6y""(x,t) - FeSy"{x,t) - J dx(y' ext (x)) 2 )5y" (x,t) 

o 

L 

+ ^y'Lt{x){ I dx x yl xt (x x )by(x x ,t)). (12) 



L 



o 
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Using the boundary conditions for the hinged end points, we find the normal modes of the rod, 5y(x, t) = 5y n (x)e luJnt . 
At the saddle point, y e xt(x) = yo(x) = 0. Using this in Eq. ([12|) . we obtain 

/2 T17T 
— sin(— at), (13) 

with n=l,2,3.... The normal mode frequencies at the saddle point are given by (J is used to denote the saddle point) 



^Lnear.n = w « \/ti 2 - — . (14) 



= TTi/ — . n = 1 is the unstable mode and it has the imaginary frequency u>% = iVLnneari where 



Q Linear — 1^*0 \j~ !• (15) 



For small amplitude vibrations around the buckled state, one has to put y ex t(x) — Ay — sin( j^x). The normal modes 
are the same as in the Eq. (??), but the normal mode frequencies are different. They are 



uj n = uiQny n 2 — 1 for n > 1, (16) 

while 



u;, =-„ r '2( j-\ ). (17! 



The Rate near first buckling instability 

The reaction rate using quantum transition state theory may be written as (lij 

c\fric oc 

jyS ric — n f 1 ^ 1 \ r L Linear \ TT j£n »-/3AgM";f!' r /io\ 



where, we have used the following expression for fi£^ c ear [111 



2tT ^Linear n CJr ■ 

n— 2 Ltnear.n 



ti L Zear = -l + xl(l) 2 < "M 



2 V V2/ u U 

and the expression for p is given by 



l 



Pi 



^ L Linear T T / Linear ,n \ J — 1 



n 



oo / | iv i \ 

J] 2 I 2sinh(ifi/3 X *) ^ 2s inh(ift/3A* ~) J 

oo / ! W \ ' 

n I 2siiih(±ft/3x„) P 2sinh(ift/9A„,j) J 



(20) 



In the above ui n and 4 represents frequencies of normal modes of the system at reactant geometry and at the 
saddle point respectively in absence of bath. In case when the system is coupled to the bath A* is the new 'bath' 
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frequency and xfi ^ s the new 'system' frequency at the saddle point, also X n j is the new 'bath' frequency and Xn is 
the new 'system' frequency at the reactant geometry. To evaluate p, we use the following two identities [1 81 ] 

oo 2 

s inh(x)^xY[(l + 1 ^) (21) 
fe=l 

and 

oo 2 

sm(x)=xl[(l-- j !L s ). (22) 
fe=i 

Following Pallak [l6| one can prove the following identity 

JV oo ( JV 

i-^Zar + a I1( A M + ^ >< II ^ + ^ I1( A 5 + ^ 

i=l n=2 V j=l 

iV 00/ JV \ 

- (-oL„ ear + e 2 + er) n ^ 2 + f2 ) x n + ^ n ^ + a ), m 

a=l n=2 \ a=l ) 

where e is an arbitrary number. Also 

TV oo I N 



(Xi 2 + ^ W^l, + e 2 ) x [] (X 2 + c 2 ) I[( A «.,- + e2 

i=l n=2 \ j=l 

/V oo / JV 

= (a;? + , 2 + e7 ) n k + a x n + e2 ) n (- 2 + e2 ) i • ( 24 ) 



We use the notation v = || . Insertion of Eq. into Eq. (J5TJ) and use of Eq. ((23} and Eq. with e = gives 

f ft (kV + xl) \ ft ft f' < tV + A ») N i 

ft ( (ft <<■ v + ^ ft ft f < fc v + X 1A \ \ r,,> 

Interchanging the order of the products and using ratio of Eq. (|23j) and Eq. (|24|) with the identification e = kv gives 
the desired result. 



n 



^ (fc 2 ^ 2 - 2 „ c + vkl) J 

n f n ) . (26) 

i3vi.i(* 3 ^+«Sn«r,»+''*7)y 



On simplification we get 



r 

Pf = - 



7+2^-^/ 7 2 +4^|„ ieai , \ / 7+2t/+^/ 7 2 +4£l 



7+2^-^/ 7 2 -4(^i) 2 \ p / 7+2^+V7 2 -4(a>l) s 



X 



7+2»>- v /7 2 -4(a J * I „ ear| „) 2 \ / 7+2^+ A /7 2 -4(^* 1TleaT , -T1 ) 2 
2t/ I I 2v 

TT _V / V L (27) 

p f 7+2^-V7 2 -4(^) 2 \ p / 7+2^+V7 2 -4(^) 2 ^ 
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BEYOND THE SECOND BUCKLING INSTABILITY 

As y^f- ~~ ¥ 2, u\ inear2 — > and the reaction rate diverges [see Eq. (18)]. So the rate expression in Eq. (18) is valid 
only if one is not too near ./— = 2. This is due to the setting in of the second buckling instability. As the rod 



is compressed, first the mode A\Jj^ sin(|-a:) becomes unstable and this is the first buckling instability and the rod 
buckles as a result of this. The length at which this occurs shall be denoted by Lf. If one supposes that the rod is 



compressed further keeping the straight rod configuration, then at a length L s , the mode A 2 y sin(-^-x) too would 
become unstable and this is the second buckling instability. What happens here is a reaction path bifurcation for the 
crossing from one buckled state to the other and is very interesting. For e > 4e c , there is only one saddle point but 
for e < 4e c , this saddle point bifurcates into two and consequently the calculation of rate near the bifurcation is a 
challenging problem. In a similar fashion one can have the third instability at a length L t etc. but these present no 
problem as far as rate calculation is concerned (see below). In order to analyze the rate near and beyond the second 
buckling instability, we proceed as follows. We assume that the Eq. (12) has solutions of the form 

y (x) =A 1 J- sm(-i) + A 2 J - sm(—x). (28) 

Using this, the elastic potential energy is given by 

Ftt a {A\ + AA 2 ) 2 Fir 2 A 2 (e-e c ) 



V[A U A 2 ] =- 



8L 4 L 2L 2 
2Fir 2 A 2 {e-Ae c ) 
L 2 ' 



(29) 



The two buckling instabilities are clearly evident from this expression - as each the coefficient of A\ or A\ changes 
sign from positive to negative. Finding the extrema of this potential leads to the following solutions for [A\, A 2 ) 

1. (0, 0): this is the straight rod configuration. Between first and second buckling (i.e. L s < L < Lf ), this is 
the saddle point. But after the second buckling, it is no longer a saddle, but it becomes a hill top. It has an energy 

Ehilltop — 0. 

2. (±^y/LL (e c — e), 0): These are the buckled states and both of them have the same energy Ef, = — ^f k {e — e c ) 2 . 

3. (0, ±i^/L_Lo(4e c — £))'■ These are the two new saddle points that arise from the bifurcation of the one that 
existed for 4e c > e. These two have the same energy Ef^ le = —^-^{e — 4e c ) 2 . At these saddle points, the rod has 
a bent (S-shaped) geometry. Beyond the second buckling instability, the barrier height is given by 

AEE:? ner = _M^££(_ 2£ + 5£c ). (30) 



Beyond the second buckling instability, away from the instability, one can do a normal mode analysis near the vicinity 
of the new saddles - there are two of them, both making identical contributions to the reaction rate. Near the buckled 
state, the normal modes have the frequencies given in the Eq. (fT6| and Eq. (fT7|) . while near the saddle, the frequencies 
are given by 

^Bent.l = ^Bent, (31) 



4e„ t ,2 = <W8(^-4), (32) 



U% Bent,n = W" 2 " 4, (33) 



for n > 2. In the above w Bent , and has an imaginary frequency with 



&Bent — VitOQ. 



(34) 
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Now the quantum rate beyond the second buckling instability can be calculated taking the saddle to be the bent 
configuration. 

Ctfric oo 

= 2 x p a (^)&=*) II (35) 

n—2 Bent,n 



where, we have used the following expression for ^i B "^t II] 

)/ric _ _ 



-2+VQ) +3^ (36) 



and the expression for p s is given by 

r f ^+ 2<y -\ 

2t/ ' ^ I 2l/ 



7+2^- A / 7 2 -4(^i) 2 \ p ( 7 +2^+V7 2 -4(^i) 2 



r / T+^-fa^^.JM p / 7+2^+ v /7 2 ^4(4 e „ t J2 \ 

x TT _V / V _ / (37) 

We have multiplied the quantum rate by a factor of 2 to account for the fact that there are two saddles of equal energy. 
It is interesting that the normal modes for this saddle retain their stability irrespective of what the compression is. 
The first mode is always unstable and other modes always stable for all values of e. Therefore, this rate expression is 
valid for all values of e < 4e c - that is even through the third buckling instability. 



The rate near second buckling instability 



Near the second buckling instability (yf t — ► 2), both oj\ inear 2 and u Bent 2 vanishes, causing the rate to diverge [see 

Eq. (18) and Eq. (35)]. The cure for this divergence is to go beyond the harmonic approximations for the first two 
modes. Our discussion here follows that of Weiss jl2|. For the model described by the Hamiltonian 



H = ^Jd Xf iy 2 +V[y] + \Y,j dx 



We note that m a {x) is independent of x and may be written as m a , lo\ (x) too, as well as c Q (x). So we can write the 
above expression as 



If 2 

H=- dx/j,V +V[y} + 
o 

J L ( 



c a y(x) 



2 



x)+m Q ^Ua(x)-^^ ). (39) 
Now let 

^w = EVz^" sin(! ^F ) (40) 
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and 



^) = EVi A " sin(! ^ ) - (41) 



Then, 

H = pi^Vn +F[Ai,A 2 ,....A n ] 



+ O E i TO « ^a,n +™<*^a £a,n - — — J 



(42) 



which decouples all the modes. The Hamiltonian for the first two modes coupled with bath is given by 



1 V-r S ? 2 ( ^ C » A 1 V 

a \ sa,l ™ , ,2 / 



A \ 2 



+ Ua,2-^i } (43) 

In the Euclidean action (S^) contains contributions from the system (Sg), the reservoir (S^') and the interaction 



/3 ft 

S S E) + J dT(L^+L^) (44) 



with 

/3 ft 



+ ^(^) 2 + ^ (:?2)(A2(T) - :?2)2 

+ L 2 (A 2 )}. (45) 
In this optimized action two modes are decoupled. Also 

N 

2 



N 1 2 2 

4 B ' = E o™«&«,1 + ^/ + £ Q ,2 + b&a/) (46) 



Q = l 



4^ - E(-C6,,l^l + - CataflA* + \^_) (47) 



Now we write Euclidean action for the two modes (with bath) separately. For the first mode 

^ - ■ 2 





W 1 2 1 r 2 4 2 
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and for the second mode 



ME) [ , A ( dA 2 (T) s ' 



+ ^u 2 2 (A 2 ){A 2 {r) - A 2 f + L 2 (A 2 ) 

N 



+ £4m-(£, a +«fe 2 2 ) " ^,2^2 + ^^y}} (49) 

a— 1 " 

( E) 

The expression for 5} is equivalent to the expression for an inverted parabolic potential coupled to N harmonic 
oscillatore. So we have a system of (7V + 1) harmonic oscillators and we can calculate JV" C using the method discussed 
in the classical calculation [l(|. Now we will calculate the partition function using the 'action' defined in Eq. (|49[) . 

Following Weiss [l2| we define "influence action" s\^j(A 2 ,£), which captures the influence of the environment on the 
equilibrium properties of the open system. 

4 £) = S<$[A*€\ = 4 E) £] + S^IA^Z] (50) 
The stationary paths of action, which we denote by A 2 and by £, a ^ 2 , obey the Euclidean classical equations of motion 

-f- dV(A 2 ) " 

a — l 

m o^a,2 ~ m a^l£, a ,2 + c aM = (51) 

We choose for convenience to periodically continue the paths £, a ,2(i~), A 2 (t) outside the range < r < j3Ti by writing 
them as Fourier series 

n—oo 

fa,2(T)= £ U2,ne iUnT 

72— — CO 
71— CO 

A 2 (t)= J2 A 2 ^ T (52) 
n=— oo 

where S, a ,2,n — £a,2,-n; ^2,n = ana1 v n = jfir * s a bosonic Matsubara frequency. Substituting Eq. (f52|) into Eq. 

O and Eq. ijlS]). we obtain 

iV n— oo 

a— 1 n— — oo 



(l£|£a,2,nf 



2 

T^2.n 



(53) 



Next, we decompose £ a ,2,n into classical term £ a 2 n and a deviation y a ,2,n describing quantum fluctuations, 

-^2,« + y Q ,2,n- (54) 



£a,2,n — £ Q ,2,n + Va,2,n 
Ca 



In the second form, we have used the solution of the oscillator mode £ a ,2,n following from Eq. (fSTj) . Since £a,2( T ) 
is a stationary point of action, the term linear in the deviation is eliminated and we find an expression in which the 
quadratic forms of y and A 2 are decoupled, 

41 [A2,l + y] = Sjf] [y] + ^ [M (55) 
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Quantum Rate at T=0.01K 
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FIG. 3: Plot of logarithm of rate of crossing from one buckled state to the other using quantum transition state theory 
without friction (solid line) and with friction (dashed line), for a Si rod of dimensions Lo — 500 A,w = 20 A, d = 10 A, 
7 = 1 x 10 n sec _1 and T = 0.01 K. For this rod the first three instabilities occur at L f -L = -0.1646 A, L s -L = -0.65971 
and 7 = 1 x 10 11 sec -1 . 



N n—oc 



a—1 n— — oo 

N fm 

= E / dT ^(yl.2 + u 2 a yl. 2 ), 

a=l J n 1 



(56) 



4> 2 ] = E^E 



2m a 

a—l n— — oo 



(57) 



12 



In the Fourier series representation, the influence action (Eq. I|57p) takes the complete form 

n— oo 

4f/J^] = ^ £ (c«i^,„i 2 ), 



(58) 



n— — oo 



where 



i N r 2 

L \ " " 

^ ^1 2TOq 



2 1 
/17T /3?i 



2 + <) 



Assuming Ohmic friction 12] 



Cn =Vnl- 



(59) 



(60) 



We want to calculate the ratio of partition functions for the second mode at the saddle and at the reactant geometry. 
As all other modes are harmonic, there contributions to the rate is calculated using the same procedure as discussed 
in the section 3.1. In our calculation we use the 'action' for the second mode at the saddle point, given in the Eq. 
P9"| . Expanding the kinetic term into its Fourier components, we get 



s { 2 E) = fipn 



+ u 2 2 {A 2 ) + i/ n7 ) \A 2 , n \ 2 + L 2 (A 2 ) 



(61) 



So the partition function for the second mode at the saddle point is given by 



Ql 



I*9(A 2 ) 

dA 2 e ft% 4tt 



r [^(^+7-V7 2 -4^(A 2 ))] 



x r[g(^ + 7 + y/<f-4u*(A 2 ))] x Q,„ 



th- 



(62) 



We have already determined £! 2 , oj 2 (A 2 ) and L 2 (A 2 ) variationally in the absence of reservoir. For the second mode 
at the reactant geometry, the expression for the partition function is given below 



Q-2 



dA 2 e 2kT J3h Air ., 
r [^(S + 7+AA 2 -4c 2 2 )]xQ ba 



Hi ■ 



l 4tt /37I 

In the regime where L > L S1 the rate may be calculated using (transition state is assumed to be straight rod) 



rifric oo 



PJ" — Q ( )( 

quantum^ rf V c\ ) V 



i n ea r ^2 



(63) 



(64) 



n=3 Linear,n 

In the regime where L < L s , the rate may be calculated using (transition state is assumed to be bent rod) 



R 



fric 

quantum,! 



ofric oo 
l-K It Bent W2 



(65) 



n— 3 ^Bent,n 
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Quantum Rate at T=30K 
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FIG. 4: Plot of logarithm of rate of crossing from one buckled state to the other using quantum transition state theory without 
friction (solid line) and with friction (dashed line), for a Si rod of dimensions Lo = 500 A,w = 20 A, d — 10 A and T = 30 K. 
For this rod the first three instabilities occur at Lf — Lo = —0.1646 A, L s — Lo = —0.6597 A and 7 = 1 x 10 sec -1 . 

RESULTS 

We now use the above formulae to calculate the rate of passage from one buckled state to the other, for Si rod of 
dimensions L = 500 A, w — 20 A, d— 10 A, the case considered by Carr et al. 0. The Young's modulus and density 
of Si are Q = 130 GPa and p = 5000 kg.m~ 3 . The first and second buckling instabilities occur at lengths which 
we denote as Lj and L s respectively and their values are given by Lf — Lq = —0.1646 A, L s — Lq = —0.6597 A. In 
the regime L > L s , the saddle point has straight rod configuration and in the regime L < L Sl the saddle point has 
bent configuration. At T — 0.01 K we do not show the rate calculated near the second buckling instability - as it is 
extremely low (Fig. The equations for rates have product over all the infinite modes, as would be normal for a 
continuum model. This however is not realistic, and this can be easily modified to account for discreteness of these 
rod. We do this, by restricting the product to contain just the same number of normal modes as the perpendicular 
degrees of freedom for the discrete model. We have thus taken contributions from first 2128 normal modes of the 
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rod. Fig. [4] shows the quantum rate against compression, made at a temperature of 30 K. Our calculation including 
friction shows that friction lowers rate of conversion from one buckled state to the other. 



SUMMARY AND CONCLUSIONS 



In this paper we have discussed the buckling of a nano rod under compression and rate of it's conversion from 
one buckled state to the other using quantum version of multidimensional transition state theory. As compression 
increases, the rod buckles (first buckling instability), and has now two stable states. From one stable state it can go 
over to the other by thermal fluctuations or quantum tunneling. Using a continuum approach, we have calculated the 
rate of conversion from one state to the other using system plus reservoir model. The saddle point for the change from 
one state to the other is the straight rod configuration. The rate expression, however, diverges at the second buckling 
instability. At this point, the straight rod configuration, which was a saddle till then, becomes hill top and two new 
saddles are generated. The new saddles have bent configurations and as the rod goes through further instabilities, 
they remain stable and the rate calculated according to the harmonic approximation around the saddle point remains 
finite. However, this rate too, diverges near the second buckling instability. In the quantum transition state theory 
calculation we have calculated centroid partition function for the second mode, and derived expressions that are well 
behaved through the second buckling instability. Using these expressions, we have calculated the rates for nano-rods 
of dimensions L = 500 A,w = 20 A, d — 10 A. 
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